Magnetization and specific heat of TbFe 3 (1303)4: Experiment and crystal-field 

calculations 
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We have studied the thermodynamic properties of single-crystalline TbFeg(B03)4. Magnetization 
measurements have been carried out as a function of magnetic field (up to 50 T) and temperature 
up to 350 K with the magnetic field both parallel and perpendicular to the trigonal c-axis of the 
crystal. The specific heat has been measured in the temperature range 2-300 K with a magnetic field 
up to 9 T applied parallel to the c-axis. The data indicate a structural phase transition at 192 K 
and antiferromagnetic spin ordering at Tn ~ 40 K. A Schottky anomaly is present in the specific 
heat data around 20 K, arising due to two low-lying energy levels of the Tb 3+ ions being split by f-d 
coupling. Below Tn magnetic fields parallel to the c-axis drive a spin-flop phase transition, which 
is associated with a large magnetization jump. The highly anisotropic character of the magnetic 
susceptibility is ascribed mainly to the Ising-like behavior of the Tb 3+ ions in the trigonal crystal 
field. We describe our results in the framework of an unified approach which is based on mean-field 
approximation and crystal-field calculations. 

PACS numbers: 75.30.-m, 75.40.Cx, 71.70.-d 



I. INTRODUCTION 

The family of ferroborates with the general formula 
RFe3(B03)4 (R is a rare-earth ion or Y) attracts con- 
siderable attention since its members exhibit a wide 
variety of phase transitions. They possess magnetic 
structures which change as a function of temperature, 
external magnetic field and substitutions in the rare- 
earth subsystemjii^i 3 -^ It was recently discovered that 
RFe 3 (B03)4, (R = Gd, Nd) exhibits multiferroic fea- 
tures, i.e. the coexistence of elastic, magnetic and electric 
order parameters The various ordering phenomena 
and their interaction cause anomalies in the dielectric 
permeability, electric polarization and magnetostriction, 
both spontaneous and field-induced. Our study aims to 
elucidate this interplay for the ferroboratc TbFc3 (603)4, 
in which both the rare-earth ions and the Fe-ions form a 
magnetic subsystem. 

The crystal structure of TbFc3(B0 3 )4 at high tem- 
peratures is trigonal and belongs to the space group D\ 
(R32)i&2 The main elements of the structure are spiral 
chains of edge-sharing FeOg octahedra running along the 
c-axis. The rare-earth ions are coordinated by triangu- 
lar R06 prisms which are isolated from each other since 
they are separated by regular BO3 triangles and have no 
common oxygen ions. Both the BO3 triangles and ROg 
prisms connect three FeOe chains. There arc no direct 



Fe-O-Fe links within the same c-plane and the shortest 
interchain exchange paths are therefore given by Fe-O- 
R-O-Fe and Fc-O-B-O-Fe, while the Fe-O-Fe exchange 
paths within the chains are much shorter (see Fig. [1} . In 
the case of GdFe 3 (B0 3 )4, the room temperature struc- 
ture was found to be transformed to the Df (P3i2i) space 
group at T s = 156Ki.£ Below T s , two non-equivalent Fe 
sites are present, which give rise to both a stretched and 
compressed modification of the FeOg spiral chains. 

In GdFe3(B03)4, the presence of different magnetic 
sublattices causes several magnetic transitions below Tg, 
which are associated with antiferromagnetic spin order- 
ing of the Fe-sublattice (Tn ~ 37 K) and a spin reori- 
entation at Tsr« 9K4 ' 10 i 1:L i 12 In addition, strong mag- 
netostrictive and magnetoelectric effects below Tsr are 
associated with a spin reorientation which occurs in ex- 
ternal magnetic fields with B\\ c-axis 

In this paper, we report on the thermodynamic proper- 
ties of TbFe 3 (B0 3 ) 4 . Our data on TbFe 3 (B0 3 ) 4 show a 
spin-flop phase transition, which is driven by magnetic 
fields B\\c. This transition is associated with a large 
jump of the magnetization. The highly anisotropic mag- 
netic susceptibility is ascribed mainly to the Ising-like 
behavior of the Tb 3+ ions in the trigonal crystal field. 
Our experimental data are analyzed in the context of 
the unified approach which is based on mean-field ap- 
proximation and crystal-field calculations. The spin-flop 
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FIG. 1: The crystal structure of RFe 3 (B0 3 ) 4 . The spiral 
chains formed by FeOe octahedra along c-axis (the shortest 
Fe-O-Fe exchange) are represented by the chains of shadowed 
areas. 

transition is apparently accompanied by magnctoelastic 
effects. Magnctoclcctric effects might therefore be associ- 
ated with the spin-flop transition, leading TbFe 3 (603)4 
to be a candidate for a multiferroic compound. 



II. EXPERIMENT 

The single crystals of TbFe 3 (603)4 were grown using 
a Bi 2 Mo 3 Oi2-based fluxJ^ The seeds were obtained by 
spontaneous nucleation from the same flux. Single crys- 
tals were green in color and had a good optical quality. 
The orientation of the crystals was performed by x-ray 
diffraction. All magnetic measurements were performed 
with the external magnetic field either parallel or perpen- 
dicular to the c-axis of the crystal. The AC-susceptibility 
at 1000 Hz and DC-susceptibility were measured in mag- 
netic fields of 0.001 T and 0.1 T, respectively, in the tem- 
perature range 1.8-350 K using a Quantum Design Physi- 
cal Property Measurement System (PPMS) . The magne- 
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FIG. 2: Specific heat in zero magnetic field (a) and mag- 
netization at B — 0.1 T (b) of TbFe 3 (B0 3 )4. The inset (c) 
highlights the Schottky anomaly of C p /T at T ~ 20 K. T s and 
Tn label the (presumably) structural phase transition and the 
onset of antiferromagnetic spin ordering in the Fe sublattice, 
respectively. 

tization M(B) was measured between 4.2 and 120 K in 
fields up to 15 T with a vibrating sample magnetome- 
ter (VSM).— The field sweep rate was approximately 
0.2T/min. In addition to this, a pulsed field magne- 
tometer was used for magnetization measurements up 
to 50 T.— Here, the total pulse duration amounted to 
0.05 s. The temperature dependence of the specific heat 
was measured in the temperature range 1.8-300 K and in 
magnetic fields up to 9 T with a PPMS calorimeter. 



III. EXPERIMENTAL RESULTS 

TbFc 3 (1303)4 exhibits two phase transitions below 
room temperature, as is illustrated by the specific heat 
and magnetization data in Fig. O By comparison with 
GdFe 3 (B03)4, the first order phase transition at T$= 
192 K is tentatively attributed to a structural symmetry 
reduction to the low symmetric P3i2i-phase.— Actually, 
the data exhibit several first order anomalies of the spe- 
cific heat in the vicinity of Tg which are probably related 
to structural changes. The origin of these anomalies can 
not be clarified by our present study which focuses on 
the magnetic properties of TbFe3(B0 3 )4. The structural 
phase transition is associated with an entropy jump of 
AS" w 1.9 J/mol K. At T N = 40 K, the onset of antifer- 
romagnetic spin order in the iron subsystem is demon- 
strated by a sharp jump in the specific heat and a kink 
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in the magnetization. We observe a specific heat jump of 
AC p « 59.4 J/mol K. This result agrees with the jump 
AC P which is predicted from the mean field theory for 
the antiferromagnetic spin ordering of the 5 = 5/2 Fe- 
sublattice>i& 

AC p = S ! S + {S {S + +l )2 R = 59 J / mo1 K ' W 

with R being the gas constant. The entire spin entropy 
of the Fe 3+ subsystem appears to develop only below 
Tn- The Fe 3+ subsystem must therefore be considered 
as a classical 3D antifcrromagnet, since short range spin 
correlations at significantly higher temperatures which 
would be present if the Fe 3+ chains were able to form 
quasi one-dimensional magnets are clearly absent. An 
analysis of the specific heat data^ indeed confirms that 
the magnetic ordering only occurs in the Fc subsystem 
while the Tb subsystem is polarized by the Fe subsystem. 

In contrast to the anomalies due to the structural and 
magnetic phase transitions at Tg and Tn, respectively, 
the Schottky anomaly in the specific heat at ~20K re- 
veals the temperature-driven population of the ground 
state of the Tb 3+ ions split in the magnetic field of the 
Fe subsystem. No spin reorientation (similar to the one 
seen in GdFe3 (1303)4 ferroborate) was observed at low 
temperatures in zero magnetic field. 
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FIG. 3: Magnetization, at T = 4.2 K, for the magnetic field 
applied to the c- and a-axis, respectively. Data obtained in a 
quasi-static field (VSM) are shown in (a) whereas (b) displays 
pulsed magnetic field data (see text). 

Depending on their direction, external magnetic fields 
can have a drastic effect on the magnetic properties of 
TbFe3 (1303)4. As shown in Fig. [3J for applied mag- 
netic fields parallel to the a-axis there is a linear field 
dependence of the magnetization at 4.2 K. The data im- 
ply Xa=0.14/is/T f.u. In contrast, applying a magnetic 
field along the easy magnetic axis, i.e. the trigonal c-axis, 



causes a sharp jump in the magnetization at -B c (4.2 K) ~ 
3.5 T. Here, the magnetic field drives a first order tran- 
sition. We demonstrate further below that the main fea- 
ture of this transition is a spin-flop of the AFM ordered 
Fe 3+ spins which is accompanied by the orientation of 
magnetic moments of the Tb 3+ subsystem along the ex- 
ternal magnetic field. The strongly discontinuous charac- 
ter of the transition is illustrated also by the pulsed field 
data in Fig.^b). In the pulsed field (dB/dt ~ 5000 T/s) 
the hysteresis at the spin-flop transition is one order of 
magnitude larger than in the quasi-static measurement 
(dB/dt ~ 10~ 3 T/s) shown in Fig. ^a), due to the re- 
laxation processes and magnetocaloric effect. At high 
magnetic fields B > B c , the data show a linear M(B) 
curve with Xe=0.13 /^s/T f.u. being slightly smaller than 

Xa- 

The magnetization data in Fig. [3] show a strong 
anisotropy. It is reasonable (and will be shown in Sec. 
[V] in a detailed analysis) to attribute this anisotropy to 
the highly anisotropic Tb 3+ magnetic moments, which 
are subject to the staggered field of the antiferromagnet- 
ically ordered Fe spins. The discussion hence must not 
only address the two subsystems of Fe and Tb moments, 
of which the Fe system divides in the two sublattices be- 
low Tn, but also the Tb subsystem must be discussed 
in terms of two sublattices since it is polarized by the 
indirect exchange interaction with the iron subsystem. 

The field dependences of the magnetization at 4.2 K 
shown in Fig. [31 reveals two anisotropy features: (1) A 
sharp increase of the magnetic moment ~9^b evolves 
at the critical field Bq. (2) The linear contribution to 
M(B\\c) only exists above Bq. If the extra moment (1) 
is neglected, the linear part of M(B\\c) is a straight line 
through the origin which is very similar to M(B\\a). Such 
behaviour (2) is typical for an uniaxial antifcrromagnet 
magnetized along the easy axis. 

Upon heating (see Fig. [4j, the spin-flop transition 
shifts to higher temperatures and vanishes above Tn, but 
it remains clearly visible in the entire antifcrromagncti- 
cally spin ordered phase. However, the magnetization 
jump reduces in value and becomes less sharp as it is 
shifted towards higher fields. For fields larger than Be 
as well as for T >Tn the magnetization curves become 
nonlinear. 

Below Tn and for fields B\\c <Bq, the staggered field of 
the Fe spins polarizes the terbium subsystem due to f-d 
interactions yielding two Tb sublattices with the mag- 
netic moments oriented oppositely and along the c-axis. 
For B\\c the effective field acting on the terbium sublat- 
tice which magnetic moment opposes the applied field is 
decreased as the applied field increases and this magnetic 
moment tends to diminish. This is a factor determining 
the character of the magnetization at 4.2 K prior to the 
phase transition, since at low temperatures the longitu- 
dinal susceptibility of the iron subsystem is very small. 
At Be the staggered field acting on the Tb 3+ magnetic 
moments vanishes as the magnetic moments of Fc 3+ ions 
become oriented nearly perpendicular to them. As a re- 
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FIG. 4: (Color online) Magnetization for the magnetic field 
applied to the c-axis at various fixed temperatures. 

suit, the magnetic moments of Tb 3+ ions align with the 
external magnetic field reaching the saturation value of 
about 9 \ib at low temperatures. 

In the spin-flop phase both terbium sublattices have 
magnetic moments directed along the field B\\c and the 
magnetic moments of iron sublattices bend towards the 
field direction. This part of the magnetization curve al- 
lows us to estimate the transverse susceptibility of the 
iron subsystem: \± e ~ Xc = 0.13 /is/T f.u. The ob- 
servation of Xa > Xc suggests a small contribution from 
the Tb subsystem for B\\a. At higher temperatures, as 
shown in Fig. 4, the magnetization curves gradually be- 
come less sharp. In the initial collinear phase the iron 
subsystem begins to contribute to the magnetization as 
longitudinal susceptibility grows with increasing temper- 
ature, thus stabilizing the initial phase. 

The magnetic phase diagram, constructed from this 
data and from specific heat measurements in magnetic 
fields, is presented in Fig. [5] The difference in phase 
boundaries shown by open and full symbols accounts 
for the hysteretic phenomena and finite slope of M(B) 
curves at the spin-flop transition. This diagram shows 
that magnetic fields parallel to the c-axis slightly sup- 
press the antifcrromagnctic spin ordering temperature of 
the Fe sublattice. In contrast, the critical field of the 
spin-flop transition increases upon heating. 



IV. THEORY 

The magnetic properties of terbium ferroborate are 
governed by both magnetic subsystems interacting with 
each other. As already mentioned above, the crystal 
structure of rare-earth ferroborate*!^^ suggests that 



FIG. 5: (Color online) (a) Specific heat anomaly at Tn for 
different magnetic fields and (b) Magnetic phase diagram for 
B\\c from magnetization and specific heat data. Bq labels 
the critical field of the spin-flop transition. Full (open) data 
markers refer to increasing (decreasing) quasi-static fields, i.e. 
dB/dt ~ 1(T 3 T/s (see the text). 



the main feature of their magnetic structure is formed 
by the helical iron chains along the c-axis. However, our 
experimental data provide no evidence for quasi-lD mag- 
netic properties, which indicates significant interactions 
between the chains. Moreover, the experimental magne- 
tization curves and phase diagram have the form which is 
typical for 3D uniaxial antiferromagnets. This fact pro- 
vides reason enough to consider the magnetic structure of 
the Fe subsystem as that of iniaxial aitiferromagnet with 
two sublattices magnetic moments of which are equal in 
the absence of field. The interactions hierarchy in the sys- 
tem is such that it stabilizes the orientations of Fe and Tb 
magnetic moments along the trigonal axis. Indeed, the 
strongest interaction in the compound is the crystal field 
(the splitting of the ground multiplet comprises about 
400-500 cm -1 , the energy interval between the ground 
quasi-doublet and the nearest excited states is of the or- 
der of 200 cm" 1 ), and the non-Kramers Tb 3+ ion in the 
trigonal crystal field acquires a huge anisotropy (g c ~ 18 , 
g a ~ ) and becomes an Ising ion with the preferred 
orientation of magnetic moment along the trigonal axis. 
Interaction comparable in value with the crystal field is 
absent in the system, therefore there is no source of any 
other orientation of Tb magnetic moments. And so there 
are grounds to suggest a collinear arrangement of mag- 
netic moments in TbFc3 (603)4 . We note however that 
only neutron diffraction studies on an HB-enriched com- 
pound can unambiguously resolve the magnetic structure 
under discussion. 

In the following section we model the magnetic proper- 
ties of TbFc3(BC>3)4 using molecular field theory Within 
this approximation the total Hamiltonian, which com- 
prises the Hamiltonians of the iron and the terbium sub- 
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systems as well as the one considering the Tb-Fe interac- 
tion, can be expressed as a sum of single-particle Hamil- 
tonians. In the presence of an external magnetic field B, 
the effective Hamiltonians of the Tb/Fe ion of the i-th 
sublattice (i=l and 2) can be written as: 

Hi(Tb) = W CF - sjMbJJB + B mi (Tb)], (2) 

Hi(Fe) = -g sf x B Si[B + B TOi (Fe)], (3) 

where gj is the Lande factor and Jj is the angular mo- 
mentum operator of the rare-earth ion; 175 = 2 is the 
g-factor and Si is the spin momentum operator of the 
Fe ion; B m j (Tb/Fe) are the molecular fields acting on 
the Tb or Fe ion in the z-th sublattice. The crystal-field 
Hamiltonian TLcf is governed by the symmetry of the 
local environment of the rare-earth ion. 

As explained in Scc. lIIIi from the analogy to the struc- 
tural phase transition at 156 K in GdFe3 (1303)4 we sup- 
pose that the specific heat anomaly at T§ = 192 K re- 
flects a crystal symmetry reduction from the trigonal 
space group R32 to the trigonal group P3i2i. This im- 
plies that the local symmetry of the Tb 3+ ion is reduced 
from D 3 at T > T s to C 2 at T < T s . Nevertheless 
(as will be discussed below) a simpler Hamiltonian of D$ 
symmetry is sufficient for describing the low-temperature 
magnetic properties of TbFe3(B03)4, since the available 
crystal-field parameters of isostructural compounds give 
rise to two closely-spaced (with a splitting less than 2-3 
cm -1 ) singlets as a low- lying state of the ground mul- 
tiplct. These singlets are responsible for the Ising-likc 
behavior of this non-Kramers ion in both the paramag- 
netic and the ordered phases. 

In order to describe the thermodynamic properties 
of rare-earth compounds we only need to consider the 
ground multiplet. In terms of the equivalent operators 
O™ the crystal field Hamiltonian TLcf reads as follow- 
ing: 

Hcf = ajBlOl + (3j{BlOl + BlOl) (4) 
+7./(i? 6 O 6 + B 6 3 O 6 3 +i? 6 6 O 6 6 ), 

where S™ are the crystal-field parameters for the D3 
symmetry, aj, (3j and jj are the Stevens coeffi- 
cients. The former are unknown for the Tb 3+ ion in 
TbFe 3 (1303)4. However, our results depend only weakly 
on the actual crystal-field parameters since the magnetic 
properties of TbFe3(B03)4 are essentially controlled by 
the magnetic behavior of the iron subsystem and by the 
Ising character of the non-Kramers Tb 3+ ion in the crys- 
tal field of the trigonal symmetry. We hence could not 
determine the crystal-field parameters from the exper- 
imental data. Instead, we have used the crystal-field 
parameters of the isostructural compounds, in particu- 
lar the rare-earth aluminoborates (see e.g. Ref. [l7l ) and 
refined these using the data for crystal-field splitting of 
Nd 3+ in NdFe 3 (B0 3 ) 4 ^ 



We have written the expressions for the molecular 
fields acting on the Tb and Fe ions from the i-th sublat- 
tice on the basis of the interactions hierarchy and mag- 
netic structure of TbFe3 (603)4. A detailed analysis of 
the interactions which was performed for the low tem- 
perature phase of GdFc3(B03)4 in Ref. [l^ also seems to 
be quite applicable to TbFc3 (603)4. According to this 
analysis, a rare-earth ion does not interact with iron ions 
from the same c-plane but interacts antiferromagnetically 
with iron ions from the neighboring planes, it is clearly 
shown in Fig. 7b of Ref.12. The interaction between rare- 
earth ions can be neglected because the distances be- 
tween them are larger than 6 A and no reasonably su- 
perexchange path could be introduced. Consequently, 
none of the ferroborates demonstrates an intrinsic order- 
ing in the rare-earth subsystem. The molecular fields 
B m i(Tb) and B m i(Fe) can therefore be written as 

B mi (Tb) =X fd Mi, (5) 

B roi (Fe) = XMj + X fd m l , j = 1 or 2, j ^ i, (6) 

where Xfd < and A < are the molecular constants 
of the Tb-Fe and Fe-Fe antiferromagnetic interactions. 
The magnetic moments and of the i-th iron and 
terbium sublattices are defined as 

M l = 3gsHB(Si) and m, = gj(i B (Ji)- (7) 

The thermodynamic potential of the system (per for- 
mula unit) has the following form: 

1 2 

$(T,£) = -[-fcsT^ln^Tb) (8) 

i=l 

I 2 

+5 I>WJi)B mi (Tb) 

i=i 

2 

-3fc B T^lnZi(Fe) 

i=i 

1 2 

i=l 

with the partition functions 

Z i {Th,Fe) = ^M ~ En f b 1 ! Fe \ 0) 

where E n i(Tb/ Fe) are the eigenvalues of corresponding 
Hamiltonians (Eqs.2 and 3). To find the magnetic mo- 
ments of the Tb and Fe subsystems Eqs.7 in the mean- 
field approximation one has to solve the self-consistent 
problem of deducing their values and orientations on the 
basis of Hamiltonians Eqs. [2] and [3] while considering 
the condition of the minimum of the thermodynamic po- 
tential (Eq. [S]) for a given temperature and field. The 
right part of equation for Mi is the relevant Brillouin 
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function, as it should be in the case of an equidistant 
spectrum which is typical for the Fe 3+ ion with an or- 
bital singlet as a ground state (S-ion). Substituting the 
magnetic moments for particular phases into the thermo- 
dynamic potential (Eq. [8]), we obtain the energies of the 
phases and find the critical fields for the phase transitions 
from the condition of their equality. The magnetization 
of the compound (per formula unit) is given by: 



i 2 



(10) 



Both magnetic subsystems contribute to the initial 
magnetic susceptibility of terbium fcrroboratc: 



XP + XP, Ml or 



(11) 



In the paramagnetic range, where the interaction be- 
tween iron and terbium subsystems can be neglected, the 
magnetic susceptibility of the Tb subsystem can be calcu- 
lated with the help of the well-known Van Vleck formula 
on the basis of the crystal-field Hamiltonian Eq. 2J For 
the ordered phase at T < TV as well as for the para- 
magnetic phase, the initial magnetic susceptibilities of 
the compound can easily be found from the initial linear 
parts of the magnetization curve M(B) calculated for the 
field along and perpendicular to the trigonal axis. 

The contribution of the Tb 3+ subsystem to the heat ca- 
pacity of the TbFe3 (1303)4 compound can be calculated 
with the help of the usual expression (per rare-earth ion, 
i.e. per formula unit): 



C T b = fcj 



(E 2 ) - {Ef 
(k B T) 2 



(12) 



The thermal averages were calculated for the spectrum of 
Tb 3+ ion formed by the crystal field and by interactions 
with the iron subsystem and an external magnetic field. 



COMPARISON OF EXPERIMENTAL DATA 
AND THEORETICAL CALCULATIONS 




FIG. 6: (Color online) Experimental (symbols) and calculated 
(lines) magnetization curves of TbFe3(BOa)4 for various tem- 
peratures with B || c. 



A. Magnetization 

In order to carry out a quantitative analysis of the mag- 
netization curves for TbFe3 (603)4 in accordance with 
the approach presented in the previous section, the mag- 
netizations in the collinear and flop phases were calcu- 
lated in the mean-field approximation and the critical 
fields of the phase transitions Be were found from the 
equality of their thermodynamic potentials. In Fig. [5] 
the experimental and calculated magnetization curves are 
displayed for several temperatures. 

The analysis of the magnetization curve at 4.2 K 
provides a possibility to determine some parameters of 
the compound. The f-d interaction parameter A/d = 
— 0.253 T//is was calculated by fitting the initial part of 



M(B) where the contribution of Fe subsystem is neg- 
ligibly small because of its small longitudinal suscep- 
tibility at low temperatures. This value of Xfd is in 
good accordance with the experimental data for the low- 
temperature splitting A of low-lying singlets of the Tb 3+ 
ion in TbFe3 (1303)4 which is around 32cm~ 1 J£ Reason- 
able crystal-field parameters for Tb 3+ in TbFe3 (603)4 
yield the g-tensor component along the trigonal axis 
in the range from 17.5 to 17.8, resulting in a value of 
A(T = OK) = [i B g z B jd (T = OK) = (i B g z \\ fd \M(T = 
OK) « 31cm- 1 . 

The value of the exchange Fe-Fe parameter, responsi- 
ble for the transverse susceptibility, Xi = — 3T//zb, has 
been obtained from the slope of the M(B) dependence 
measured at T = 4.2 K in the flop phase. However, if 
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FIG. 7: The B — T phase diagram for the spin-flop transition 
in TbFe 3 (B0 3 )4 at B || c (cf. Fig.gp ). The symbols and the 
line represent the experimental data and the calculated phase 
boundary, respectively; the difference in experimental phase 
boundaries reflects a finite slope of M(B) curve caused by a 
sample shape and hysteresis phenomena. 



calculated magnetization curves for all temperatures as 
A2 = —2T//-LB- The B — T phase diagram for the spin- 
flop transition (see Fig. [7]) demonstrates a good agree- 
ment between the calculated and experimental data. We 
emphasize that quantitative description of the magneti- 
zation isotherms of TbFe3 (1303)4 in the mean-field ap- 
proximation requires an introduction of two exchange 
parameters characterizing Fe-Fc antiferromagnetic inter- 
action, Ai and A2. The Ai enters Eq. [8] for the energy 
and determines the critical field of the spin-flop transition 
Be and the transverse susceptibility of the iron subsys- 
tem in the flop phase. The A2 determines the magnitude 
of iron magnetic moments Mi for specific temperature 
and field as well as the Neel temperature. The necessity 
of introducing two exchange parameters is a consequence 
of the chain structure of the ferroborate when consid- 
ered in the mean-field approximation. The mtrachain 
exchange interactions are described by Ai, whereas the 
mierchain exchange interactions (which are responsible 
for the three-dimensional magnetic order) are described 
by A2 . Our experimental data and the results of our the- 
oretical analysis indicate that these two quantities are of 
the same order of magnitude. 



the exchange parameter for the collinear phase is taken 
equal to this value then the calculated phase transition 
critical field exceeds the experimentally observed value 
and a similar relation between experimental and calcu- 
lated critical fields takes place for all temperatures. On 
the other hand, if one takes the exchange parameter in 
the collinear phase to be ~ 1% less than that in the flop 
phase, the energy of the collinear phase increases slightly, 
the critical field decreases and its calculated value is in a 
good agreement with the experimental one (see Fig. [S]). 
The difference between the values of the exchange pa- 
rameter A for two phases could be explained by the pres- 
ence of the magnetoelastic effects accompanying the field- 
induced phase transition. In support of this hypothesis, 
it is knownS that, in GdFe 3 (B0 3 ) 4 with the S-ion Gd 3+ , 
the field-induced phase transition in the low-temperature 
phase (T < 10 K) results in a magnetostriction jump 
of the order of 2 • 1CT 5 . In TbFe 3 (B0 3 ) 4 the highly 
anisotropic terbium ion should give a much more essential 
contribution to the magnetostriction at the field-induced 
spin- flop transition than the S-ion Gd 3+ . Indeed, calcu- 
lations of the multipolc moments Q nm = a ra (0™) of the 
Tb 3+ ion in TbFe 3 (B0 3 )4 (which provide an information 
on the rare-earth contribution to the magnetostriction^) 
have shown that their changes at the spin-flop phase tran- 
sition are very significant. 

To analyze the temperature evolution of the magne- 
tization we have to take into account the exchange be- 
tween the Fe chains which is responsible for the three- 
dimensional ordering. This value of A enters the Bril- 
louin function in Eq.7, the well-known explicit form of 
which is not given for the sake of brevity. It was cho- 
sen from the best agreement between experimental and 



B. AC magnetic susceptibility 

As has already been discussed above, the strong 
anisotropy of the magnetic susceptibility is caused by 
the contribution of the Tb 3+ ion with its Ising charac- 
ter of magnetic behavior. The x(T) dependencies have 
been calculated using parameters of TbFe 3 (B0 3 )4 de- 
rived from the analysis of the magnetization isotherms. 
However, as one can see in Fig. [8] the experimentally ob- 
tained and calculated XJ-(T') curves differ substantially. 
The non-zero susceptibility x\\ at very low temperatures 
and a growth of x± upon increasing the temperature up 
to Tn could be ascribed to a misorientation of the sam- 
ple by a few degrees with respect to the crystallographic 
axis. In fact, the inclusion of a misorientation by 3° 
with respect to the basal plane leads to a significantly 
better agreement between the calculated and experimen- 
tal data. The variation of the crystal-field parameters 
within reasonable limits does not result in any significant 
changes in the description of the x(^) dependencies; the 
agreement between experimental and calculated curves 
remains roughly the same. 

A distinctive convexity of the experimental X||(^) 
curve in the temperature range 30-35 K is associated with 
a Schottky-type anomaly, i.e. with the change of occu- 
pation of two low-lying terbium singlets. This anomaly 
can also be observed in the calculated curves, especially 
in the susceptibility of the Tb subsystem. The calcu- 
lated Neel temperature is higher than the corresponding 
experimental value but this is an usual feature of the 
mean-field theory. 
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FIG. 8: (Color online) Temperature dependencies of the ini- 
tial magnetic susceptibility in TbFe3(BOa)4 along (x||) and 
perpendicular (x_i_) to the trigonal axis, (a) Symbols are ex- 
perimental points; solid thick lines are calculated curves for 
X|| and xx; thick dashed line represents calculation of \± f° r 
the misorientation angle with regard to the basal plane in 3° . 
(b) Thin solid and dot-dashed lines are calculated curves for 
iron and terbium subsystems, respectively. 



C. Heat capacity 

In the following, the specific heat and the effect of the 
magnetic field will be analyzed in terms of our approxi- 
mation. The experimental specific heat data in magnetic 
fields B — and 3T along the c-axis are presented in 
C/T vs. T coordinates in Fig. [9^,. As one can see, the 
magnetic field shifts the Schottky anomaly to lower tem- 
peratures, specifically for C(T) from ~ 20K at B = to 
~ 5 K at B = 3 T. The anomaly at around 5 K is partic- 
ularly evident since in this temperature range both the 
phonon and magnetic contributions are small. 

The calculated temperature dependencies of the Tb 
contribution to the specific heat are presented in Fig. [5}d. 
The calculated Crr^T) for B = demonstrates the 
Schottky anomaly around 20 K. As has already been 
mentioned above, the anomaly is caused by the occupa- 
tion of the singlet level separated from the ground singlet 
by the temperature-dependent splitting due to the Tb-Fc 
interaction (Fig. fTOj) . 




T(K) 



FIG. 9: (a) Specific heat of TbFe3 (1303)4: open and closed 
symbols represent experimental data at B — and 3 T respec- 
tively. Thick and thin solid lines correspond to calculations 
of the Tb contributions for B = and 3T respectively, (b) 
Calculated specific heat of the Tb subsystem: dashed lines 1 
and 2 show calculated contributions for different sublattices 
at B = 3 T. 

Our analysis in Fig. [9Jd suggests that the Tb contri- 
bution to the specific heat is different for the two Tb 
sublattices if an external magnetic field is applied along 
the c-axis. For the terbium sublatticc with magnetic 
moments opposed to the field direction, the Schottky 
anomaly shifts to lower temperatures. This is easy to un- 
derstand because an increase of the magnetic field leads 
to a decrease of the energy gap between two singlets (see 
Fig. rrrj|) . It means that an external field demagnetizes 
this terbium sublattice. For the other terbium sublat- 
tice, however, the energy gap between the two singlets 
increases, the increased external field magnetizes the sub- 
lattice and the Schottky anomaly moves to higher tem- 
peratures. We emphasize the agreement between the ex- 
perimental and numerical data which confirms that the 
description of Tb-Fc coupling in TbFc3(BOs)4 provided 
in this work is appropriate. 



VI. CONCLUSIONS 

Our measurements of the specific heat and the mag- 
netization and the theoretical description of the experi- 
mental data provide a detailed understanding of the thcr- 
modynamical properties of TbFc3 (1303)4. Our unified 
approach describes the presence of the spin-flop transi- 
tion for magnetic fields applied parallel to the trigonal 
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FIG. 10: Temperature dependencies of the splitting of the 
ground quasi-doublet of the Tb 3+ ion in TbFe3(BOa)4 at B = 
and 3 T for sublattices with magnetic moments (1) opposite 
to the field direction and (2) along the field. 

c-axis and the related magnetization jump below Tjj. Mi- 
croscopically, the highly anisotropic magnetization data 
have been attributed to the Ising character of the Tb 3+ 
ion in the trigonal crystal field. In addition, the presence 
of a Schottky anomaly in the specific heat data around 
20 K and its shift as a function of the applied magnetic 
field have been explained as the Tb contribution to spe- 
cific heat in terms of two low-lying energy levels of the 
Tb 3+ ions being split by f-d coupling. 

In particular, our approximation considers the ex- 
change interaction of the Fe-subsystem as well as the 
coupling between Fe- and Tb-subsystcms. The / — d 



coupling parameter has been determined and its value 
is in good agreement with available spectroscopic data. 
The antiferromagnetic exchange interaction in the iron 
subsystem has been described with the help of two pa- 
rameters. The parameter Ai is responsible for bend- 
ing the iron magnetic moments in the flop phase, which 
occurs against the intrachain exchange, and determines 
the spin-flop transition field; the parameter A2 is con- 
nected with the interchain exchange interaction respon- 
sible for the three-dimensional ordering, it controls the 
Neel temperature and magnitudes of the iron magnetic 
moments for specific temperatures and fields. The rela- 
tion I Ai| > IA2I is naturally fulfilled because the Fe-O-Fe 
exchange paths within the chains are much shorter than 
the interchain exchange paths. The appearance of two 
exchange parameters is a special feature of consideration 
of magnetic properties of the compound with chain struc- 
ture in the framework of the mean-field approximation, 
which is known to be of limited validity for description 
of magnetic lattices of lowered dimensionality. A sub- 
stantial influence of the magnetoelastic effects which ac- 
company the field-induced phase transition on the value 
of the molecular parameter responsible for bending the 
iron magnetic moments in the flop phase has been found. 
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